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ABSTRACT 

We have made sensitive milliarcsecond-resolution radio images of the nearby 
merger galaxy Arp 299 at four epochs spread over 18 months between 2003 and 
2005. The combined data revealed a total of 30 point sources in the two pri- 
mary merger nuclei. Twenty- five of these are found in the northeastern nucleus 
(component "A"=IC 694) over a region ~ 100 pc in diameter, while 5 are in the 
southwestern nucleus (component "B1"=NGC 3690) within a region ~ 30 pc in 
size. These objects are interpreted as young supernovae and supernova remnants; 
the ratio of the source counts in nuclei A and Bl is approximately equal to the 
ratio of their predicted supernova rates. An approximate luminosity function 
has been derived for nucleus A, and indicates that it might contain as many 
as 500-1000 compact radio sources more powerful than Gas A; the integrated 
flux density of these sources would be about 20% of the total flux density seen at 
lower resolution. A new supernova occurred in nucleus Bl in the first half of 2005, 
having a peak radio power at least 2,000 times the present power of Cassiopeia 
A. This supernova is located within 0.4 pc (projected distance) of an apparently 
older supernova remnant, making it very likely that this indicates the presence of 
a massive super star cluster within nucleus Bl. Comparison of the typical radio 
flux densities of our compact radio sources to the observed X-ray luminosities of 
nuclei A and Bl indicates that it is possible that one radio source in each nucleus 
actually could be associated with an active galactic nucleus rather than being a 
supernova remnant. 

Subject headings: galaxies: starburst — galaxies: evolution — galaxies: individ- 
ual (Arp 299, NGC 3690, IC 694) — supernovae: general — radio continuum: 
galaxies 
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Introduction 



It now appears likely that the most massive galaxies in the universe formed from mul- 
tiple mergers of smaller galaxies; those mergers are apparent in the disturbed morpholo- 
gies of many galaxies image d in deep optical and i nfrared fields such as t he Hubble Deep 



Field s (IFerguson et al.l I2OOOI ) , the GOODS survev (IDickinson et al.l l2003l : iGiavalisco et al. 



20041 ). and the COSMOS survey JScoville et aP boOTh . Nearby and present-day analogi 



of these early galaxies often are found in Luminous and Ultraluminous Infrared Galaxies, 
many of which are undergoing major mergers at the present epoc h. The nearest majo r 
merger galaxy is generally regarded as the 'Antennae,' NGC 4038/9 (jwhitmore et al. 1999 ). 
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1). The starbursts in these 



mergers and other galaxies obey the so-call ed "radio-infra red relation" , whereby their radio 
and far- infrared powers are well-correlated (ICondorull992l ): the infrared emission is thought 
to be powered ultimately by the photons from hot young stars, while the radio emission 
typically is powered by synchrotron radiation from particles accelerated in supernovae (SNe) 
and supernova remnants (SNRs). 

In recent years, compact radio sources have been used as one window providing in- 
sights iiito_Jh^_jia|ure_o^ In the mo st nearby starbursts, such as 
M82 JMuxIow et all ^99^: iFenech et al1l2008h and NGC 253 Juivestad fc Antonuccil [l997l : 
Lenc &: Tingayll2006l ). interferometers reveal a mix of thermal and nonthermal radio sources 
representing H II regions and fairly young SNRs. However, in the somewhat more dis- 
tant merger galaxies, typically at distances of a few tens of megaparsecs, use of the Very 
Long Baseline Interferometry (VLBI) technique is required to separate individual compact 
sources; the interferometer sensitivity on continent-scale baselines is insufficient to detect H II 
regions, so only non-thermal emission from young S Ne and SNRs is detected. For example, a 



series of observati ons of the merger galaxy Arp 220 (IRovilos et al.ll2003l : iLonsdale et al.ll2006 



Parra et al.l 120071 ) has detected at least 49 objects in the two merger nuclei, provided the 
radio spectra of many of these objects, and detected the explosion of several new supernovae. 

The subject of our investigation is the relatively nearby (41 Mpc) m erger galaxy Arp 29 9 
(Mrk 171), which was noted as a starburst galaxy over 25 years ago by lGehrz et al.l (119831 ). 
The two primary nuclei of the merger are the northeastern component A (often referred to as 
IC 694) and the southwestern component B (often referred to as NGC 3690). Although we fo- 
cus primarily on the starburst in this paper, one or both nuclei also have bee n inferred to con 



tain active galaxy components based on X-ray imaging and spectroscopy (IDella Ceca et al. 



- 3 - 



2002 



2006|) 



Zezas et all2003l : lBallo et al.ll2004j ). optical integral-field spectroscopy ( Garcia-Marm et al 
, and the presence of nuclear H2O masers (IHenkel et al.ll2005l : iTarchi et al.l 120071 ). 



There is substantial CO emission associated with the Arp 299 starburst, and a relatively 
high ratio of HCN/CO at the two p rimary nuclei, imp l ying the presence o f dense and warm 
molecular clouds in the two nuclei ( Aalto et al.l 119971 ). iBaan fc Haschickl ( 19901 ) found OH 
megamaser emission from nucleus A, and inferred that this emission comes from a clumpy 
molecular medium located 100-600 pc in front of the conti nuum radio source (s) b eing a mpli- 



fied. Near-infrared imaging and modeling of Arp 299 by lAlonso-Herrero et al.l (I2OOOI ) have 
shown that nucleus A accounts for 50% of the near-infrared emission of the merger and 
has a star- formation and supernova rate approximately six times those of nucleus B. The 
latter nucleus actually is resolved int o a dominant component Bl, and a weaker component 



B2 lo cated an arcsecond to the NW (ILai et al.lll999l : lAlonso-Herrero et al.ll2000l : iNeff et al. 



2004 ). In order to provide context to the reader, we reproduce two previously published 
figures here. Figure 1 is an 814 nm Hubble Space Telescope image of Arp 299, while Figure 2 
is a high-resoluti on Very La r ge Ar ray (VLA) image of the merger at radio wavelengths, both 
reproduced from lNeff et al.l (120041 ). 



Arp 299 has hosted at least five optic al supernovae in the last 20 years, the last of 
which was the type II supernova SN 2005U (jMattila et al.ll2005l : iModjaz et al.ll2005l ). These 
optical supernovae all occurred outside the primary merger nuclei; those nuclei are likely to 
be the most frequent sources of supernovae in Arp 299, but their extreme obscuration makes 
finding optical supernovae highly unlikely. Even a two-epoch NICMOS near-infr ared search 
failed to reveal any new supernovae in the two merger nuclei (ICresci et al.l 120071 ). However, 
high-resolution VLBI imaging of Arp 299 rev ealed five compa ct radio sources in nucleus A, 
including at least one very young supernova (INeff et al.ll2004j ). 



This pap er reports VLBI observations of Arp 299 with considerably higher sensitivity 
than those of iNeff et al.l (|2004j ). aimed at probing deeper into the radio luminosity function 
and searching for the onset of radio emission from young supernovae. 



2. VLBI Observations 



The merger system Arp 299 was observed on four occasions using an interferometer 
including the Very Long Baseline Arr ay (VLBA) augme nted by the Robert C. Byrd Green 
Bank Telescope (GBT). The VLBA jNapier et al.lll993h is an interferometer consisting of 



10 identical 25m antennas distributed from Hawaii to the US Virgin Islands, while the GBT 
is a 100m x 110m single aperture in West Virginia, which was added to the VLBA in order 
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to approximately double the overall sensitivity of the array. Observations were carried out 
under program code BU027, and took place at approximately six-month intervals between 
late 2003 and mid-2005. 

Each of the four observing epochs consisted of almost identical 10-hr sessions. At each 
epoch, blocks of 60-110 minutes were alternated between sky frequencies of 2.2715 GHz 
(hereafter 2.3 GHz) and 8.4215 GHz (hereafter 8.4 GHz). There were four blocks at 2.3 GHz 
and three blocks at 8.4 GHz; 2.3-GHz observing blocks were placed at the beginning and end 
of each 10-hr session in order to minimize the impact of atmospheric water vapor fluctuations 
in low-elevation observing at 8.4 GHz. At each band, four 8-MHz frequency channels were 
used at both right and left circular polarizations. All were sampled at the Nyquist rate with 
two bits per sample, providing a total data rate of 256 Mbit s~^. These data were recorded on 
a mix of instrumentation tapes and hard-disk modules which were shipped to Socorro, New 
Mexico for processing on the VLBA correlator. Only the parallel polarization hands were 
processed, in order to keep the data rate within the correlator output limits; each 8-MHz 
frequency band was further subdivided into 32 or 64 spectral channels in order to preserve 
the capability for wide-field imaging. 

The GET was typically missing for 10%-20% of the time in each observing session, due 
to software and pointing problems as well as 5-10 minute periods to change between 2.3 and 
8.4 GHz. The Mauna Kea VLBA antenna did not observe in the first and third epochs due 
to bad weather. The Kitt Peak VLBA antenna had 25% of its 2.3 GHz data removed in all 
observations, due to bad interference in one of the 8-MHz bandpasses. Additional data losses 
were typically no more than 5%-10% of the data from 1-2 antennas for a given observing 
run and frequency band. 



The observations were carried out in phase- referencing mode ( iBeasley &: Conwaylll995l ). 
Observations of 2.5 minutes in length on Arp 299 were interleaved with observations of 
approximately 45 seconds on the calibration source J1128+5925, located 0.86° from Arp 299. 
Both the antennas and the correlator were targeted midway between the strongest centimeter 
radio component, nucleus A of Arp 299, and the possible backgrou nd source D, which were 



previously detected in the VLBI observations of iNeff et al.l (120041 ): this pointing preserved 
signal/noise in the relatively small beam of the GBT as well as minimizing the effects of 
bandwidth smearing. Occasional 2.5-minute observations of Jl 127+5650, located 2.58° from 
J1128+5925, were made in order to provide a check of the phase-referencing quality. Table 1 
is a summary of the observations, including the observation dates and the total integration 
time on Arp 299 at each frequency band. 



a: 



1 data calibra tion was carried out in NRAO's Astronomical Image Processing System, 



AlPS (lGreisenll2003l ). Standard corrections included those for Faraday effects in the Earth's 
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ionosphere, application of the most accurate Earth orientation parameters, gain corrections 
to remove the effects of non-optimal sampler levels, and calibration of the feed rotation of 
the individual antennas. Initial coarse solutions for the interferometer delays were provided 
by fitting one or two scans of a strong calibrator (usually J0927+3902=4C 39.25), while final 
delay and phase calibration (or "fringe fitting" ) was performed by least-squares fitting on each 
45-second scan of Jl 128-1-5925. The resulting delay and phase solutions were interpolated in 
time to provide the final calibration of Arp 299 and the check source. 

Particular care was taken in setting the amplitude calibration scale for the GBT. VLBA 
antennas provided useful system temperature data at roughly 1-minute intervals, and their 
calibration is judged to be accurate at the 5% level. However, the a priori and real-time 
calibrations of the GBT were less accurate; good values for GBT are important because it 
is the most sensitive telescope in the observing array. Nominal values for GBT sensitivity 
were used to weight the GBT during the fringe-fitting process. Then the amplitude scale was 
adjusted so that the flux densities of the amplitude-check sources J0854-I-2006 (OJ 287) and 
J1310-I-3220 varied smoothly with baseline length; in this adjustment, care was taken to make 
the amplitude on interferometer baselines to the GBT similar to those involving the relatively 
nearby VLBA antenna in Hancock, New Hampshire. This is a subjective process because 
source structure is convolved with the fringe pattern for each individual baseline, so the final 
calibration was determined "by eye" using the source flux-density plots. Adjustments of the 
GBT sensitivity relative to the nominal values were as high as 25%, but more typically in 
the 10% range. 



3. Imaging, Source Identification, and Flux-Density Measurements 

3.1. Imaging 



For each data set, images were produced using the standard deconvolution and "Clean" 
algorithms in AIPS. Initially, six fields were imaged i n each dat a set, centered on the six 
strongest and most compact VLA sources detected by iNeff et al.l (120041 ). These six regions 
were the two primary nuclei of the galaxy merger, A and Bl; source B2 to the NW of Bl; 
the two possible additional merger nuclei, C and C'; and the possible background source 
D. Each field consisted of a l'.'024-square image centered on the position of the appropriate 
compact VLA source, with at least four pixels per resolution element. The data sets were 
Cleaned to minimum fiux-density levels of 3-4 times the predicted noise for each image, in 
order to assure that sidelobes due to imperfect sampling were removed completely. 

Careful source identification procedures, discussed in Section [3^ were carried out on all 
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six fields. The only significant detections in individual epochs or in combined data from all 
four epochs were found in the fields A, Bl, and D. Therefore, to avoid generating spurious 
sources by Cleaning noise spikes, we produced the final interferometer images by targeting 
only the three fields with detected sources; this improved the noise characteristics of the final 
images by ~ 10%. The noise levels from the final three-field images are cited in Table 1, 
and are typically 10%-20% above the theoretical values under ideal conditions. The effective 
noise levels were somewhat higher because of imperfect phase referencing, as discussed in 
Section 13. 3[ 

After making the initial images, we experimented with self-calibration in order to refine 
the atmospheric calibration and hence improve the final images of Arp 299. However, be- 
cause there is very little correlated signal present, adequate signal/noise required averaging 
times considerably longer than the atmospheric coherence time, so self-calibration was not 
effective. Thus, we depended instead on the data with only the phase calibration relative to 
J1128+5925. 

The purposes of the four epochs of observation were twofold-both to achieve a time base- 
line long enough to see source flux density changes, including the appearance of new sources, 
and to achieve sufficient integrat ion time to go significantly deeper than the VLBA+GBT 



observations by iNeff et al.l (120041 ). For the latter purpose, it is necessary to combine all the 
data from four epochs. The traditional ways to do this are to (1) combine the calibrated data 
from all four epochs in order to make a single set of images, or (2) average the four-epoch 
images in the image plane. The former method may be susceptible to errors in the case 
of significant source variability, while the latter method does not easily retain the relative 
weighting of the different baselines or of observations with different noise and dynamic range 
levels. We chose the first method, combining the calibrated data sets and then producing 
three-field images of the combined data as for the individual epochs. Figures 3 and 4 show 
the final 2.3 GHz images of the regions of nuclei A and Bl containing all the detected com- 
pact radio sources. Images of the same fields at 8.4 GHz are not shown because of the smaller 
number of sources detected in that band. From these figures, we note that the compact radio 
sources in nucleus A extend over a region ~ 100 pc in diameter, while the compact sources 
in nucleus Bl are all found within a region ~ 30 pc in diameter. 



3.2. Source Identification 

The final source identification within the fields of A, Bl, and D proceeded in an iterative 
process. First, all individual peaks above five times the rms noise were identified in the 2.3- 
GHz images, and peaks above six times the rms noise were identified at 8.4 GHz, both 



- 7- 



at the individual epochs and in the combined data sets. (The higher noise threshold at 
8.4 GHz is caused by the fact that the 8.4 GHz images have 16 times as many pixels as 
the 2.3 GHz images.) At some individual epochs, there were apparently spurious sources 
above the cutoff ffux densities near the strongest (~ 1 mJy) sources, almost certainly due to 
the imperfect phase-referencing. These spurious sources were easily identified because they 
occurred at different locations at the different epochs, and largely disappeared when data 
sets were combined. 

Locations of all sources above the cutoffs in the images from the combined data sets were 
inspected at each individual epoch. Sources exceeding the nominal cutoff in the combined 
images were considered real if individual peaks above 2.5a also could be identified in at least 
two of the four epochs; alternatively, if a peak above 7a that was not an obvious image defect 
could be found in a single epoch, a source was considered real and possibly variable. In fact, 
using these criteria for the individual epochs, all sources above the cutoffs in the combined 
data sets were found to be real, and we judge the source identification to be complete at the 
5(7 (2.3 GHz) and 6a (8.4 GHz) levels for the combined data sets. 

3.3. Source Flux Densities and Positions 

For weak source detections, some sources can appear shghtly resolved because of noise 
plateaus lying near the actual sources. In our images, no sources were obviously resolved 
when this effect was taken into account; given unresolved sources, the flux densities may be 
measured by a variety of methods. We chose to make quadratic fits to the inner few pixels of 
each source in order to determine the peak fiux density, and used this value as the total fiux 
density for the apparently unresolved sources. An alternative method of fitting beam-sized 
gaussians to each source makes assumptions about fiux-density distributions that may not 
be warranted, particularly in the situation where imperfect phase-referencing scatters flux 
non-randomly in the image plane. 

As stated previously, the phase-referencing process corrects imperfectly for the tropo- 
sphere and ionosphere. Thus, there is a net coherence loss in the Arp 299 images, and an 
apparent reduction in correlated flux densities. We measured this reduction by using the 
phase-reference check source, J1127-I-5650 (hereafter J1127). J1127 was imaged using the 
initial phase cahbration apphed from J1128-I-5925, then self-calibrated and re-imaged. The 
increase in peak flux density by this procedure was a measure of the loss of coherence due to 
the phase referencing errors. The reduction in fiux density for Arp 299 at each epoch then 
was estimated by multiplying the coherence loss of J 1127 by the ratio of separations from the 
calibrator, (0.86°/2.58°), since the coherence loss should depend approximately linearly on 



- 8 - 



the separation between reference and target source (iBeasley &: Conwaylll995l ). Table 2 gives 
the flux-density correction factor for Arp 299 inferred by this method for each epoch, the 
average correction factor over the four epochs, and the "real" noise levels achieved for the 
combined data sets after the imperfect phase- referencing is taken into account. Inspection 
of Table 2 indicates that the 2.3 GHz coherence losses were more consistent from epoch to 
epoch, while the 8.4 GHz coherence losses were much worse during the northern summer; 
this supports the inference that the 2.3-GHz losses were dominated by an ionospheric cali- 
bration error whose magnitude is relatively independent of season, while the 8.4-GHz losses 
were dominated by the troposphere, which is more variable during the summer. Final flux 
densities for the Arp 299 sources were obtained by multiplying the fitted flux densities by 
the amplitude correction factor at each epoch, or by the average correction factor for the 
combined data sets. 

Errors in fiux-density estimates have three different causes. First, there is an overall 
scale error, dominated by the error in calibrating the GET, which we estimate to be 10% 
(see Section [2]). Second, there is an error in estimating the coherence loss, including the 
assumptions that this loss depends linearly on target/calibrator separation, is otherwise 
independent of direction, and has the same temporal structure for Arp 299 and J1127 (the 
latter source was observed only every hour or two). We take this error to be 50% of the 
inferred coherence loss; i.e., for an apparent coherence loss of 15%, the la error in the loss 
factor is assumed to be 7.5%. (In reality, the error in coherence loss is surely non-gaussian, 
but there is no definitive reason to adopt any other form.) Third, there is the error caused 
by the limited signal/noise ratio of the data, which is as high as 20% for a 5a source and as 
low as 3% for a 1-mJy source at a single epoch. We add these three errors in quadrature to 
derive the final errors at each epoch. For the final images using the combined data sets, we 
assumed that the fiux-density scale errors and coherence-loss errors were uncorrelated from 
epoch to epoch, implying that they are reduced by a factor of two for four epochs. 

Source position measurements are unaffected by amplitude scale errors. At 2.3 GHz, 
epoch-to-epoch position shifts should be dominated by signal/noise considerations, which 
limit the positional accuracy to approximately the beam size divided by the signal/noise 
ratio; this may be as large as 1 mas for a 5a-6a source. At 8.4 GHz, position shifts between 
epochs sometimes are larger than the 0.2-0.3 mas that this simple assumption would indicate, 
most likely due to the imperfect phase calibration. We must als o accou nt for the uncertainty 



in the position of J1128+5925; errors are given by iFey et al.l (120041 ) as 0.16 mas in right 



ascension and 0.27 mas in declination. In fact, the current VLBA calibrator lisl0 was used 



( http://ww'w. vlba.nrao.edu/astro/calib/vlbaCalib.txt ) 
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it quotes 50% larger errors, and uses updated coordinates that differ from lFey et al.l (120041 ) by 
about 0.05 mas in each coordinate. Given that the signal/noise at 2.3 GHz hmits position 
determination to 1 mas accuracy for the weakest sources, and that the phase-referencing 
may hmit 8.4 GHz positions at a similar level, we estimate position errors of 1.0 mas in each 
coordinate, and quote absolute source positions to an accuracy of 0.1 milliseconds of time in 
right ascension and 1.0 mas in declination. 

Table 3 gives the final positions and flux densities of all our confirmed sources at 2.3 GHz, 
while Table 4 gives a similar list for 8.4 GHz. In these tables, source names are given by 
convention as position offsets from aj2ooo = ll'^28™' and 5j2ooo = 58°33'; for sources in 
common between both frequency bands, the source name derived from the 8.4 GHz position 
is used. The source properties and their implications are discussed further in section HI 



4. Discussion of Compact Radio Sources 
4.1. Source Detection Summary 



Inspection of Tables 3 and 4 shows that there are significant milliarcsecond radio sources 
detected in both of the primary nuclei of the galaxy mergers, A and Bl; as stated in Sec- 
tion l3.H there were no detections in B2, C, o r C^ S ource D is detected at all epochs, but with 
no other radio detections nearby. iNeff et al.l (120041 ) previously reported five VLBI detections 
in A, as well as detection of D, but none in Bl. The much deeper observations described in 
this paper revealed 19 sources at 2.3 GHz in A, as well as 13 sources at 8.4 GHz; since only 
7 objects were detected at both frequency bands, there are now a total of 25 distinct VLBI 
radio sources in A. In addition, there were four "steady" detections in Bl, plus a variable 
source that only appeared at the fourth epoch. Since no sources are resolved at either 2.3 or 
8.4 GHz, we infer that they generally have sizes smaller than one milliarcsecond, or 0.2 pc. 
The most likely interpretation, by virtue of analogy with similar sources in other merger or 
starburst galaxies such as Arp 220, M82, and NGC 253, is that the sources are all young, 
compact SNe or SNRs. In fact, the ratio of 5:1 in source coun ts between A and Bl is ver y 
close to the ratio of 5.9:1 in their estimated supernova rates (jAlonso-Herrero et al. 2000), 
indicating that the stellar mass functions and supernova evolution may be similar between 
the two nuclei. In an extre me merger such as Arp 299 , where much of the starburst has 
occurred in the last 10 Myr (jAlonso-Herrero et al.ll2000l ). it is likely that virtually all of the 
radio sources have resulted from Type II supernovae. 
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4.2. Source Variability 



In Figure 5, we plot the four-epoch 2.3 GHz data and 8.4 GHz data for the strongest 
source at each band, one of intermediate strength, and one of the weaker sources. This figure 
includes data for only a few sources, but most of the detected sources show no significant 
variability. Since the error bars at individual epochs are typically in the 10%-15% range, 
the upper limits on variability over 18 months are only in the vicinity of 20% for the non- 
variable objects. These limits are consistent with the much smaller variability and hmits 



seen in most radio sources in NGC 253, M82, and Arp 220 (jUlvestad fc Antonuccil 11994 
I997I : iKronberg et~aDl200ol : IParra et allboOTl : IPenech et allboosf h 



We find only two sources in our program that display significant variability, both at 
8.4 GHz. The first is the strongest 8.4 GHz source shown in the right-hand panel of Figure 5, 
designated as 33.621+46.71 in Table 4; this source in nucleus A appe ars to have a syste matic 
fiux density decline over four epochs. It was initially reported by iNeff et al.l (120041 ) as a 
very young radio supernova, and found to have an 8.4 GHz fiux density of 3.2 mJy at 
epoch 2003.11; this higher fiux density in early 2003 confirms the monotonic decline seen in 
Figure 5. 

The other variable source at 8.4 GHz appeared only at the fourth epoch, in nucleus Bl, 
as noted in Table 4 (source 30.988+40.78). Further discussion of this source is deferred until 
Section 14. 5[ 



4.3. Source Spectra 

Figure 6 shows the distribution of two-point spectral indices (a, defined by Sy oc z/"*"") 
for all the detected sources except the new source that appeared at the fourth epoch; those 
sources with upper limits to their fiux densities at 2.3 or 8.4 GHz are shown hatched. Errors 
in the values of a range from 0.09 to 0.19 for the sources detected at both frequencies, while 
the errors on the spectral-index limits for sources detected at only one frequency range from 
0.16 to 0.21. Because only 8 of the 29 sources in this plot were detected at both frequencies, 
we have not attempted to find the true distribution of source spectra. Mathematically, 
one could use survival analysis to derive a hypothetical "true" spectral index distribution. 
However, this would rely on an implicit assumption that sources with upper limits at each 
frequency are distributed similarly to the sources detected at both frequencies, and this is 
very unlikely to be a correct assumption. 

We can say that at least seven of the sources have positive spectral indices, almost 
surely indicating optically thick synchrotron emission. (Brightness temperatures of these 
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milliarcsecond sources are above 10^ K for even the weakest sources, and thus the radio 
emission is highly unhkely to be thermal in nature.) The two variable sources, only one of 
which is included in Figure 6, have spectral indices a > +1.5 between 2.3 and 8.4 GHz. Such 
steep positive spectra are characteristic of very young radio SNe which have bec ome optically 



thin at the higher frequency, but remain optically thick at the lower frequency (jWeiler et al. 



2OO2I ). Hence we infer that the radio sources occupying the right side of Figure 6 are likely 



to be the most recent supernovae. 



4.4. Radio Luminosity Functions 

The distributions of radio luminosities at 2.3 and 8.4 GHz are shown in Figure 7, with 
the objects from nuclei A and Bl distinguished from one another. The detection thresholds 
are 1.7 x 10^^ W Hz"^ at 2.3 GHz and 1.8 x 10^^ W Hz"^ at 8.4 GHz. For reference, 
we compare these values to the galac tic SNR Cas si opeia A. Using the spectrum and flux- 



density ev olution of Gas A found by iBaars et al.l (Il977l ). as well as the 3.4 kpc distance 
derived by lReed et al.l (119951 ) . we derive Gas A powers at epoch 2005.0 of 1.6 x 10^^ W Hz~^ 



at 2.3 GHz and 6.1 x 10^^ W Hz ^ at 8.4 GHz. Thus, our respective detection thresholds 
for milliarcsecond radio sources in Arp 299 are ~ 11 and ~ 29 times the Gas A power at 2.3 
and 8.4 GHz. 

It is of interest to compute the total number of radio sources that might be expected 
above the Gas A power in Arp 299, as this may provide some insights into the time evolution 
of the sources. We consider only nucleus A, since the number of radio sources in nucleus Bl 
is too small to derive useful statistical results. One can use the distribution of observed radio 
powers shown in Figure 7 to fit a luminosity function, and then integrate that luminosity 
function down to the power of Gas A. However, the small-number statistics give rather large 
errors on any such fit. For a rough estimate, we assume a power-law form for the luminosity 
function, and take bin widths of 0.2 dex in radio power. At log(P) = 19.75 (in W Hz~^), we 
ameliorate small-number statistics slightly by averaging this bin with the two adjacent bins, 
and find 2 sources per bin at 8.4 GHz and 2.33 sources per bin at 2.3 GHz. A two-point fit 
between these average values at log(P) = 19.75 and the values at log(P) = 19.35 gives the 
following results, where P19 = P/(10^^ W Hz~^): 

dN{Pig)/d\nP,, ~ 80Pi9-'-' (1) 

at 2.3 GHz, and 

rfAr(Pi9)/rflnPi9 ~ 30Pi9-'-' (2) 

at 8.4 GHz. Integrating these estimated luminosity functions down to the luminosity of Gas A 
at the respective frequencies, we find that nucleus A of Arp 299 would contain approximately 
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890 sources stronger than Cas A at 2.3 GHz and 740 sources stronger than Cas A at 8.4 GHz. 
The error bars on these numbers are roughly a factor of two, based on relatively small-number 
statistics, and also ignore the possibility that the shapes of the luminosity functions might 
change at lower luminosities. Indeed, it is possible that most of the objects detected in 
our observations are fairly young radio SNe, whose emission still is largely dominated by 
interaction with their own mass-loss shells, rathe r than SNRs dorninate d by interaction with 
an external medium (jChevalier fc Franssonll200ll : IChevalier et al.ll2004l ): any transition from 
SNe to SNRs most likely would cause the luminosity function to chan ge shape. Extensiv e 
discussion of the SNe/SNR evolution in the merger Arp 220 is given by iParra et al. 
and we refer the interested reader to that paper for details that also may apply to Arp 299. 

The summed flux densities of the observed young supernovae and supernova remnants 
in nucleus A (Tables |3] and H]) amount to 5.3 mJy at 2.3 GHz and 4.2 mJy at 8.4 GHz. 
Multiplying the above luminosity functions by power and integrating down to the power of 
Cas A, we find respective total flux densiti es of 19 and 14 mJ y at 2.3 and 8.4 GHz, in compact 
radio sources more powerful than Cas A. iNeff et al.l (120041 ) gave respective flux densities of 
101 and 77 mJy at 4.9 and 8.4 GHz in the sub-arcsecond component of nucleus A. Thus, 
even if the estimated luminosity functions continue unbroken to the power of Cas A, the 
total flux density in young SNe and SNRs more powerful than Cas A is only ~15-20% of 
the total radio flux density in nucleus A, consistent with all observations. 

We can consider the implications of this result by comparing it with the expected su- 
pernova evoluti on in our owri Galaxy. It has been ~ 330 yr since the supernova event that 
created Cas A (jHugheslll980l ). so one can examine the (admittedly naive) hypothesis that 



all supernovae in Arp 299 evolve exactly like Ca s A. For a supernova rate of 0.65 yr ^ in 



nucleus A of Arp 299 (lAlonso-Herrero et al.ll2000l ). there should now be ~ 200 SNRs younger 
than Cas A in that nucleus, while integration of the estimated luminosity functions implies 
that Arp299-A may contain as many as ~ 800 SNRs above the Cas A power. This is weak 
evidence that the SNRs in Arp 299 stay stronger for a longer period of time, as one might 
expect if they explode in a dense medium and expand more slowly than in our own Galaxy. 
By comparison, the nearby starburst galaxies M82 and NGC 253 have ~ 1 0-20 radio SNRs 
more powerful than Cas A (jUlvestad fc Antonuccilll997l : iFenech et al.ll2008l ) in galaxies with 
supernova rates of ~ 0.1 yr~^, which would scale to ~ 60-120 objects for supernova rates of 
~ 0.6 yr^^. This implies that the SNR evolution may be more rapid, or the SNRs intrinsi- 
cally less lu r ainou s, in these weaker starbursts. Alternatively, as suggested for Arp 220 by 
Parra et al.l (120071 ). it may be that the initial mass function in the more extreme starburst 
galaxies is more top-heav y. However, the evidence for the existence of such top-heavy mass 
functions is inconclusive (lElmegreenl |2005| ) . 
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Given the large uncertainties in the luminosity function, a significantly deeper VLBI 
observation of Arp 299, which should be possible within 2-3 yr, would be required to test 
whether the young SNe and SNRs in Arp 299 are more radio-powerful or longer- lasting than 
those in less active starbursts. A turnover in the luminosity function that would be caused by 
a top-heavy initial mass function, as well as additional young supernovae, might be revealed 
by such an observation. 



4.5. The Two Newest Supernovae 



As mentioned in Section we have detected two compact variable sources at 8.4 GHz, 
which we identify as y oung supernovae . First, we consider the source 33.621-1-46.71, previ- 
ously denoted as AO by lNeff et al.l (120041 ). This object had an 8.4 GHz flux density of 3.2 mJy 
at epoch 2003.11, and fell to 1.1 mJy at epoch 2005.55. Fitting an exponential falloff of the 
form S{t) = S'o(to) (t/to)"^ since epoch 2003.11, we flnd that the source flux density at 

8.4 GHz was declining as However, since the source was not detected at 2.3 GHz at 
any of our four epochs, we conclude that it still remained optically thick at that frequency 

2.5 yr after becoming optically thin at 8.4 GHz. The slow rise time at 2.3 GHz is consistent 
with a Type II supernova, but not a Type Ib/c supernova (see Table 3 of Weiler et al. 2002). 

The other variable source was detected only at our last epoch, 2005.55, in Nucleus Bl; 
it had an 8.4-GHz flux density of 1.17 ± 0.28 mJy at that time. We have no later VLBI 
observations to study its evolution, but we did search for all relevant observations in the 
VLA data archive. We found two A conf iguration observa tions that were taken subsequent 
to the many epochs of data analyzed by et al.l (120041 ). These observations were made 
under program code AC749, on 2004NOV02 (2004.84) and 2006APR15 (2006.29), and have 
resolution of 0'.'2. We have extracted both data sets from the VLA archive and analyzed 
them. Historically, the 8.4 GHz flux densities of nucleus Bl ranged from 5.9 to 7.2 mJy in 
A conflg uration data ob tained between 1990 and 2002, while nucleus A ranged from 70 to 
75 mJy (iNeff et al.ll2004l ). Thus, a new radio source with a flux density of a few millijansky 
or less would not be noticeable within nucleus A, but might be detectable in nucleus Bl. 
Table 5 gives the results of our flux-density measurements of nuclei A and Bl from the 
archival VLA data. We flnd that the flux densities of both nuclei were consistent with their 
historical values at epoch 2004.84, whereas the 8.4-GHz flux density of Bl had approximately 
doubled in 2006.29. 



Our best estimate for the change in the flux density of Bl is an increase of approxi- 
mately 6.5 mJy be tween 2004.84 and 2006.29. Given the estimated supernova rate of only 
~ 0.11 yr~^ in Bl ( lAlonso-Herrero et al.ll2000l ). it seems reasonable to attribute this entire 
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increase to the single new source we have detected in our VLBI observations. Given that 
assumption, Figure 8 shows an 8.4 GHz hght curve for the Arp 299 radio supernova, using 
the combination of archival VLA data and VLBI observations. The supernova in Arp 299 
reached a peak power of at least 1.2 x 10^^ W Hz~^ at 8.4 GHz, about 2,000 times the Gas A 
power. This is comparable t o the p eak radio powers of the luminous Type Iln supernovae 
summarized bv IWeiler et al. ( 20021). It is fa r more powerful than the late-time powers of 
Type Ibc supernovae (jSoderberg et al.ll2006f), and nearly 10 times more powerful than the 
recently discovered radio supernova in M82 (IBrunthaler et al.l 120091 ). 

We note the estimated supernova rates of 0.65 yr^^ for nucleus A an d 0.11 yr~^ for Bl 
jAlonso-Herrero et al.l[200o[ ). Our 8.4-GHz VLBI observations and those of iNeff et all (j^OOJ) 
spanned a total of ~ 2.5 yr; detection of one new supernova in each of A and Bl is consistent 
with those values, assuming that most of the young supernovae actually are detectable radio 
supernovae. 



4.6. Supernovae as Signposts for Super Star Clusters 

The new supernova in Bl appeared within two milliarcseconds, or 0.4 pc (projected), of 
a steady 2.3 GHz source with a flux density of 0.25 mJy; this 2.3 GHz source did not vary 
when the new source appeared. The four steady 2.3 GHz sources all lie within a rectangle 
about 100 by 160 mas, with a total area of 1.6 x 10^ mas^. One can test the hypothesis that 
a new supernova occurring with equal probability at any location in this rectangle would 
appear very close to an existing source by chance. For a completely random location of 
the new supernova, the a priori probability that it would occur within 5 mas of an existing 
source is 2.0%, and the probability of occurrence within 2 mas of an existing source is 0.3%. 
Thus the likelihood is that the new supernova is physically associated in some way with the 
existing radio source. If the two radio sources actually are SNe/SNRs separated by 0.4 pc, 
their supernova shells would overlap in no more than 100 yr even for expansion speeds of 
only 2,000 km s'^ 

The most logical physical association is that both the older SNR and the new supernova 
are located in the same su per star cluster (SSG) in nucleus Bl. Such SSGs have typical radii 



of 4 pc in the Antennae (jWhitmore et al.lll999l ). correspon ding to 20 mas at Arp 299. In 



fact, the SSGs found nearby, such as R136 in 30 Doradus (IHunter et al.lll995l ). are highly 
centrally condensed, so massive stars and supernovae may be most likely to coexist over a 
much smaller region. Since the compact VLBI sources detected in Bl are otherwise sepa- 
rated by typical distances of 10-15 pc, it is tempting to hypothesize that each is a beacon 
identifying the location of an individual SSG. We note that a number of SSGs dominated 
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by thermal emission previously have been detected and imaged by the VL A (ITurner et al. 



2OOOI : iBeck et al.l I2OO2I : Ijohnson fc Kobulnickvl l2003l : iTurner fc Beckl I2OOJ). However, the 



SSCs dominated by thermal emission apparently are younger than 3 Myr, whereas the pos- 
sible SSCs indicated by our VLBl observations must be older than 3 Myr in order to harbor 
SNe and SNRs. 

Statis t ical i nferences based on single objects are, of course, quite perilous. However, 
Neff et al.l (12004 ) also noted that the new 8.4-GHz source in nucleus A appeared within 3 pc 
(projected) of a previously detected 2.3-GHz source, and suggested that they might be in 
the same SSC. There are at least four compact radio sources in nucleus A within a projected 
separation of ~ 10 pc, making it likely that these are associated in some way. A deeper 
VLBI integration would reveal whether there are weaker SNRs physically associated with 
the same regions. 

The possible confusion of radio sources at different frequency bands provides some cause 
for anxiety in considering the spectral indices of the eight objects apparently detected at both 
radio frequencies (see Section S3]), because there is no guarantee that detections at multiple 
frequencies actually correspond to the same source. The eight multi-frequency objects do 
not vary significantly at 8.4 GHz over the course of 1.5 yr. Most appear to have optically 
thin synchrotron spectra, and thus are likely to be older SNRs radiating at both 2.3 and 
8.4 GHz rather than chance coincidences. However, object 33.630+46.79 has a fiat spectrum 
between 2.3 and 8.4 GHz, so it is possible that it is an SNR with an unusual spectrum or 
that it does, in fact, correspond to two confused sources. 



4.7. Possible Active Galactic Nuclei 



In recent years, there have been several report s that one or both of Arp 299 nuclei A and 



B (or Bl) may be active galactic nuclei (AGNs). iDella Ceca et al.l (120021) d etected a hard 
X-ray component in Arp 299; arcsecond-scale imaging by IZezas et al.l (120031 ) indicated that 
nucleus Bl has a hard X-ray spectrum and may be an AGN, while s pectr al fits to nucleus 
A indicated that it too could have an AGN component. iTarchi et al.l (120071 ) also found H2O 
maser emission in both A and Bl, supporting the possibility of dual AGNs in Arp 299. 

Since we detect numerous radio sources in each nucleus, it is not possible to determine 
which (if any) might correspond to the X-ray sources. However, we can at least check 
whether the radio/X-ray ratios for our typical compact radio sources might be consistent with 
the presence of AGNs. We make use of the definition Rx = vL^ib GHz)/Lx(2 — 10 keV) 
(ITerashima &: WilsonI |2003| ) . In nucleus Bl, the X-ray source has an X-ray luminosity of 
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~ 7 X 10^^ ergs in the 0.1-10 keV range (jZezas et al.ll2003l ): our typical compact radio 
source (excluding the single young supernova) has a strength of ~ 200 /xJy at 5 GHz (based on 
the 2.3 GHz detections and spectral limits), giving vL^ ~ 1 x 10'^^ ergs s~\ and log Rx ~ —3.8 
(ignoring relatively small c orrections for the X-ray spectrum). The X-ray source in nucleus A 
emits ~ 1 X 10^° ergs s~^ ( Zezas et al. 20031 ): the relatively strong 1 mJy radio source which 



Neff et all (120041 ) speculated to be a possible AGN has pL^ ~ 3 x 10^*^ ergs s"^ at 5 GHz, 



so \ogRx ~ —3.5 if this radio source were associated with the X-ray emission, or somewhat 
lower (similar to Nucleus Bl) if one of the weaker radio sources were the X-ray source. 
Thus our estimated values of log Rx ~ —4.0 to —3.5 are con s istent with the values of low- 
luminosity Seyfert galaxies observed by iTerashima &: Wilson! (120031 ). and it is possible that 
one of the compact radio sources in either A or Bl (or both) may be an AGN. However, it 
also is possible that an AGN could have \ogRx ~ —4.5, and hence be undetected in the 
present VLBI observations. The best way to determine which of the compact radio sources 
might be an AGN may be to conduct a high-sensitivity phase-referenced VLBI observation of 
the H2O masers in order to see if they could be identified with any of the compact continuum 
radio emitters. 



Neff et al.l (120041 ) hypothesized that VLBI resolution of source D could indicate whether 



it belongs to Arp 299 or is a background AGN. None of our four epochs of observations have 
resolved this radio source, so we su spect that the po ssible resolution on the scale of the beam 
size that was reporte d previously (INeff et al.l 120041 ) may have been spurious. In fact, given 
the X-ray luminosity ( Zezas et al. 20031 ) and average radio flux density of this object, it has 
l og Ry ~ —2.3, consistent with values found for some PG quasars by ITerashima &: Wilson 
J2003f ). 



5. Summary 

We have observed the nearby galaxy merger Arp 299 at four epochs and two frequencies 
with a long-baseline radio interferometer consisting of the VLBA and the GBT. The primary 
results of these observations follows: 



1. Thirty compact radio sources were detected in the two primary merger nuclei above 
limits of 10-30 times the luminosity of Gas A. The 25 detections in nucleus A span a 
diameter of ~ 100 pc, while the 5 detections in nucleus Bl span a diameter of ~ 30 pc. 



2. The ratios of the numbers of radio sourc es in nuclei A and Bl are co nsistent with the 
ratios of the supernova rates inferred by lAlonso-Herrero et al.l (120001 ) . 
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3. Most sources have variability upper limits of ~ 20% over 18.5 months, consistent with 
the interpretation that they are relatively young supernova remnants. 

4. A few objects were detected at 8.4 GHz and not at 2.3 GHz, indicating that they may 
be younger radio supernovae that are still optically thick at the lower frequency. 

5. A previously detected young supernova in nucleus A had a steadily declining 8.4 GHz 
flux density over 2.5 yr, but still had not become detectable at 2.3 GHz after the same 
interval. 

6. A new supernova was detected in nucleus Bl, with an apparent explosion date in the 
first half of 2005; this object had a peak luminosity at least 2,000 times the Gas A 
power. The new supernova occurred within 2 mas (0.4 pc projected) of a previously 
known steady radio source; it is highly likely that these two radio sources are the 
signposts of a super star cluster that formed at least 3 Myr ago in nucleus Bl. 

7. Gomparison of the radio powers of the individual sources to the X-ray luminosities of 
nuclei A and Bl indicates that it is possible that one radio source in each nucleus is 
actually an active galactic nucleus rather than a young SNR. 



The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. I thank the 
many NRAO staff that made these observations possible. I also thank Susa r i Neff for assis- 
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of the paper. 
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Fig. 1.— HST WFPC2 814 nm image of Arp 299, taken from iNeff et all t004 ]. The 
primary merger nuclei A and B, the possible smaller nuclei C and C, and the additional 
radio source D are all labelled. Reproduced by permission of the AAS. 
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Arp 299 8.4 GHz (3.6 cm) 




11 28 34.0 33.0 32.0 31.0 

Right Ascension (J2000) 



Fig. 2. — VLA 8.4 GHz image of Arp 299, taken from iNeff et al.l (120041 ). Labels are similar 
to those in Figure 1, except that nucleus B is separated into sources Bl and B2. Reproduced 
by permission of the AAS. 
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11 28 33.65 33.64 33.63 33.62 33.61 33.60 33.59 
Right Ascension (J2000) 
Grey scale flux range= (-50, 250) uJy/beam 



Fig. 3. — VLBI image of nucleus A of Arp 299 at 2.3 GHz, combining data from all four 
epochs. The intensity scale is shown at the top. 



-24- 



1 00 200 




11 28 30.995 30.990 30.985 30.980 30.975 
Right Ascension (J2000) 
Grey scale flux range= (-50, 250) uJy/beam 

Fig. 4. — VLBI image of nucleus Bl of Arp 299 at 2.3 GHz, combining data from all four 
epochs. The intensity scale is shown at the top. 
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Fig. 5. — Time histories of strengths of compact sources at three different flux-density levels, 
in nucleus A, at 2.3 GHz (left) and 8.4 GHz (right). Each flux-density curve is labelled by the 
source designation. Two sources (33.620+46.70 and 33.640+46.58) are shown in both panels. 
The fading 8.4 G Hz source, 33.621+46.71, is the apparent supernova reported previously by 
Neff et all J2004h . 
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Fig. 6. — Distribution of two-point spectral indices of compact radio sources in nucleus A. 
Hatched values represent either lower limits (a > 0) or upper limits (a < 0), while the eight 
sources detected at both frequency bands are shown un-hatched. 
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Fig. 7. — Histogram of radio luminosities of the compact radio sources in the two nuclei of 
Arp 299. The sources in nucleus A are shaded, while those in nucleus Bl are unshaded. The 
recent supernova in nucleus Bl is not included. 
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Fig. 8. — Light curve of radio supernova 30.988+40.78 in nucleus Bl at 8.4 GHz, showing 
four VLBI observations and the VLA archival observation from 2006.29. Open circles are 6a 
upper limits from our VLBI observations, and la error bars are shown on the two detections. 
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Tablc 1. VLBA+GBT Observations of Arp 299 



(1) 

Date 


(2) 

Sky Frequency 
(GHz) 


(3) 

Integration 
(min) 


(4) 

Resolution 
(milliarcsec) 


(5) 

rms noise*^ 
(/iJy beam^^) 


2003 Dec 29 


2.2715 


178 


6.80x5.37 


36 




8.4215 


180 


1.93x1.31 


27 


2004 Jun 27 


2.2715 


178 


6.29x4.69 


29 




8.4215 


181 


1.76x1.14 


27 


2005 Jan 02 


2.2715 


172 


6.37x4.97 


27 




8.4215 


173 


2.16x1.19 


31 


2005 Jul 17 


2.2715 


174 


6.49x4.85 


35 




8.4215 


176 


1.84x1.16 


30 



^The "raw" rms noise for each individual image is given in this table. 
These values should be scaled upward by the amplitude correction factors 
shown in Table 2 to estimate the "true" noise values. 



Table 2. Amplitude Correction Factors Caused by Coherence Losses 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Frequency 


2003 Dec 29 


2004 Jun 27 


2005 Jan 02 


2005 Jul 17 


Average 


Measured rms 


Effective rms 














(/ijy beam~^) 


(/ijy beam~^) 


2.3 GHz 


1.15 


1.24 


1.11 


1.19 


1.17 


18 


21 


8.4 GHz 


1.07 


1.31 


1.23 


1.43 


1.26 


15 


19 
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Table 3. Compact Sources Detected at 2.3 GHz 



Name 


a (J2000) 
11*^28"^ 


8 (J2000) 
58°33' 


Flux Density 

(/^Jy) 


Comments 


Nucleus A 










33.594+46.56 


33.5941^ 


46'.'560 


354+32 


Also at 8.4 GHz 


33.612+46.70 


33.6124' 


46'.'696 


119+22 




33.615+46.67 


33.6154^ 


46'.'667 


125±23 




33.617+46.72 


33.6173' 


46'.'724 


115±22 




33.618+46.70 


33.6175' 


46'.'695 


382+33 




33.619+46.46 


33.6194' 


46'.'463 


149±23 




33.620+46.70 


33.6200' 


46'.'699 


1090+76 


Also at 8.4 GHz 


33.621+46.60 


33.6209' 


46'.'597 


217+26 




33.622+46.66 


33.6218' 


46'.'655 


662+49 


Also at 8.4 GHz 


33 624+46 77 


33 6241' 


46'.'771 


121+23 




33.627+46.44 


33.6267' 


46'.'440 


165+24 




33.629+46.65 


33.6290' 


46'.'647 


114+22 




33.630+46.79 


33.6301' 


46'.'786 


335±31 


Also at 8.4 GHz 


33.631+46.40 


33.6305' 


46'.'402 


241+27 




33.636+46.86 


33.6361' 


46'.'863 


139+23 




33 640+46 58 


33 6404' 


46'.'581 


141+23 


Also at 8 4 GHz 


33.644+46.63 


33.6442' 


46'.'628 


126+23 




33.650+46.59 


33.6495' 


46'.'589 


211+25 


Also at 8.4 GHz 


oo.00U+4d.04 


oo.ooUi 






Also at 0.4 LrtlZ 


Nucleus Bl 










30.975+40.83 


30.9753' 


40'.'828 


237±26 




30.983+40.87 


30.9827' 


40'.'867 


294+29 




30.987+40.78 


30.9873' 


40'.'784 


245+27 


Near 8.4 GHz supernova 


30.995+40.78 


30.9948' 


40'.'776 


332+31 


Also at 8.4 GHz 


Source D 










33.011+36.55 


33.0109' 


36'.'549 


964+68 


Background? 
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Table 4. Compact Sources Detected at 8.4 GHz 



Name 


a (J2000) 


6 (J2000) 


Flux Density 


Comments 




11^28 


58 33 






l\l 1 1 /^l 1 G A 












33 594P 


46" 560 


X r?H- 1 


Alen at 9 3 CH7 


33.599+46.64 


33.5992^ 


46'.'637 


132+22 




33.620+46.70 


33.6199*^ 


46'.'699 


371+38 


Also at 2.3 GHz 


R91 -LAP, 71 


33 691 9^ 


46"707 


Xt: X U _l L jLixJ 


Tior*! iTiiTicr Tin "v^ 




33 6918^ 


46"655 


^04+49 

f n iH- 1 H- . 7 


Al«n at 9 3 GH7 


33.628+46.62 


33.6280^ 


46'/623 


147+23 




33.630+46.79 


33.6301^ 


46'/786 


366+38 


Also at 2.3 GHz 




33 6307^ 


46"690 


9i4_i_97 




33.636+46.68 


33.6356^ 


46'/677 


141+23 




33.639+46.55 


33.6391'' 


46'.'550 


129+22 




33.640+46.58 


33.6402*^ 


46'.'580 


118+22 


Also at 2.3 GHz 


33.650+46.59 


33.6495*^ 


46'.'591 


159+24 


Also at 2.3 GHz 


33.650+46.54 


33.6501*^ 


46^537 


267+30 


Also at 2.3 GHz 


Nucleus Bl 










30.988+40.78 


30.9875^ 


40'/784 


1170+280 


Fourth epoch only 


30.995+40.78 


30.9948^ 


40'/776 


153+23 


Also at 2.3 GHz 


Source D 










33.011+36.55 


33.0108^ 


36'.'549 


1810+161 


Background? 



'^The flux densities of source 33.621+46.71 at the four new observing 
epochs were as follows. 2003 Dec 29: 1861 ± 199 /iJy; 2004 Jun 27: 
1416 ± 264 /xJy; 2005 Jan 02: 1308 ± 204 /xJy; 2005 Jul 17: 1128 + 269 /^Jy. 
These also are plotted in the upper right portion of Fig. 5. 
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Table 5. VLA 8.4-GHz Flux Densities of Arp 299 



Epoch SsA (niJy) 



Nucleus A Nucleus Bl 



2004.84 72.7 ± 3.6 7.1 ± 0.4 

2006.29 61.1 ±6.1 13.6 ±1.4 



